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Abstract 



First measurements of dihadron correlations for charged particles are presented for 
central PbPb collisions at a nucleon-nucleon center-of-mass energy of 2.76 TeV over a 
broad range in relative pseudorapidity (Arj) and the full range of relative azimuthal 
angle (A^). The data were collected with the CMS detector, at the LHC. A broadening 
of the away-side {A(p ^ n) azimuthal correlation is observed at all At], as compared 
to the measurements in pp collisions. Furthermore, long-range dihadron correlations 
in At] are observed for particles with similar (p values. This phenomenon, also known 
as the "ridge", persists up to at least | A^ | =4. For particles with transverse momenta 
(px) of 2-AGeV/c, the ridge is found to be most prominent when these particles are 
correlated with particles of pj = 2-6 GeV/c, and to be much reduced when paired with 
particles of pj = 10-12 GeV/c. 
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1 Introduction 

Measurements of dihadron azimuthal correlations ||l]-[7I have provided a powerful tool to study 
the properties of the strongly interacting medium created in ultrarelativistic nuclear collisions ||8]- 
[11]. An early indication of strong jet-medium interactions at RHIC was the absence of high- 
transverse-momentum (high-px) back-to-back particle pairs in dihadron correlation measure- 
ments [1 J and the corresponding enhancement of low-px hadrons recoiling from a high-px lead- 
ing, or "trigger", particle [3J. The recent observations of the suppression of high-px charged 
hadrons [12 J and of asymmetric energies of reconstructed jets l(T3llT4il in PbPb collisions at the 
Large Hadron Collider (LHC) provide further evidence of jet quenching, suggesting a large 
energy loss for partons traversing the produced medium. 

At RHIC, extending dihadron azimuthal correlation measurements to larger relative pseudo- 
rapidities resulted in the discovery of a ridge-shaped correlation in central AuAu collisions 
between particles with small relative azimuthal angles (|A^| 0), out to very large relative 
pseudorapidity (|A?/|) f2','6|. Although the "ridge" has been qualitatively described by several 
different models [15 -26[, its origin is still not well understood. Some models attribute the ridge 
to jet-medium interactions, while others attribute it to the medium itself. The ridge has been 
observed for particles with transverse momenta from several hundred MeV/c to a few GeV/c. 
However, the character of the ridge for even higher-px particles, as well as its dependence on 
collision energy, is still poorly understood from the RHIC results [2J. Recently, a striking ridge 
structure has also been observed in very high multiplicity proton-proton (pp) collisions at a 
center-of-mass energy of 7 TeV at the LHC by the Compact Muon Solenoid (CMS) Collabora- 
tion [27J, posing new challenges to the understanding of these long-range correlations. 

This paper presents the first measurement of dihadron correlations for charged particles pro- 
duced in the most central (0-5% centrality) PbPb collisions at a nucleon-nucleon center-of-mass 
energy {^/s^) of 2.76 TeV over a large phase space. The results are presented in terms of the as- 
sociated hadron yields as a function of pseudorapidity and azimuthal angle relative to trigger 
particles in different transverse momentum intervals. Traditionally, trigger particles have been 
utilized to represent the direction of the leading hadron in a jet, and were required to have a 
higher momentum than all the other associated particles in the jet [[21 [6l. However, as shown 
in Ref. |27| , important information can also be obtained by studying the correlation of hadron 
pairs from the same px interval, which is particularly useful when addressing the properties of 
the medium itself. The current analysis employs both approaches. This measurement provides 
a unique examination of the ridge in the most central PbPb collisions at the highest energies 
reached so far in the laboratory over a wide range in transverse momentum (2-12GeV/c) and 
up to large relative pseudorapidity ( | At] \ 4), imposing further quantitative constraints on the 
possible origin of the ridge. 

Details of event readout and analysis for extracting the correlation functions are described in 
Section |2| the physics results found using the correlations are described in Section |3j and a 
summary is given in Section |4] 

2 Data and Analysis 

The analysis reported in this paper is based on PbPb collisions at ^/s^ = 2.76 TeV collected 
during the LHC heavy-ion run in November and December 2010 with the CMS detector. The 
central feature of the CMS apparatus is a superconducting solenoid of 6 m internal diameter. 
Within the field volume are the inner tracker, the crystal electromagnetic calorimeter, and the 
brass/ scintillator hadron calorimeter. Muons are measured in gas-ionization detectors embed- 
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2 Data and Analysis 



ded in the steel return yoke. In addition to the barrel and endcap detectors, CMS has extensive 
forward calorimetry. The nearly An solid-angle acceptance of the CMS detector is ideally suited 
for studies of both short- and long-range particle correlations. A detailed description of the 
CMS detector can be found in Ref. |28|. CMS uses a right-handed coordinate system, with the 
origin at the nominal interaction point, the x axis pointing to the center of the LHC, the y axis 
pointing up (perpendicular to the LHC plane), and the z axis along the counterclockwise beam 
direction. The detector subsystem primarily used for the present analysis is the inner tracker 
that reconstructs the trajectories of charged particles with pj > 100 MeV/c, covering the pseu- 
dorapidity region \rj\ < 2.5, where ?/ = — ln[tan(0/2)] and 6 is the polar angle relative to the 
beam direction. The inner tracker consists of 1440 silicon pixel and 15 148 silicon strip detector 
modules immersed in the 3.8 T axial magnetic field of the superconducting solenoid. 

The event readout of the CMS detector for PbPb collisions is triggered by coincident signals in 
forward detectors located on both sides of the nominal collision point. In particular, minimum 
bias PbPb data are recorded based on coincident signals in the beam scintillator counters (BSC, 
3.23 < I //I < 4.65) or in the steel /quartz-fiber Cherenkov forward hadron calorimeters (HP, 
2.9 < \f]\ < 5.2) from both ends of the detector. In order to suppress events due to noise, 
cosmic rays, double-firing triggers, and beam backgrounds, the minimum bias trigger used in 
this analysis is required to be in coincidence with bunches colliding in the interaction region. 
The trigger has an acceptance of (97 ±3)% for hadronic inelastic PbPb collisions [14J. 

Events are selected offline by requiring in addition at least three hits in the HP calorimeters 
at both ends of CMS, with at least 3 GeV of energy in each cluster, and the presence of a re- 
constructed primary vertex containing at least two tracks. These criteria further reduce back- 
ground from single-beam interactions (e.g., beam gas and beam halo), cosmic muons, and 
large-impact-parameter, ultra-peripheral collisions that lead to the electromagnetic breakup 
of one or both of the Pb nuclei [29]. The reconstructed primary vertex is required to be located 
within 15 cm of the nominal collision point along the beam axis and within a radius of 0.02 cm 
relative to the average vertex position in the transverse plane. 

This analysis is based on a data sample of PbPb collisions corresponding to an integrated lumi- 
nosity of approximately 3.12 fih^^ ||30ll3T1l , which contains 24.1 million minimum bias collisions 
after all event selections are applied. 

The energy released in the collisions is related to the centrality of heavy-ion interactions, i.e., 
the geometrical overlap of the incoming nuclei. In CMS, centrality is classified according to 
percentiles of the distribution of the energy deposited in the HP calorimeters. The centrality 
class used in this analysis corresponds to the 0-5% most central PbPb collisions, a total of 
1.2 million events. More details on the centrality determination can be found in Refs. |T4ll32| . 

A reconstructed track is considered as a primary-track candidate if the significance of the sepa- 
ration along the beam axis between the track and the primary vertex, dz/cr[dz), and the signifi- 
cance of the impact parameter relative to the primary vertex transverse to the beam, dxij/ a{dxy), 
are each less than 3. In order to remove tracks with potentially poorly reconstructed momen- 
tum values, the relative uncertainty of the momentum measurement, cr{pj)/ pj, is required to 
be less than 5.0%. Requiring at least 12 hits on each track helps to reject misidentified tracks. 
Systematic uncertainties related to the track selections have been evaluated as discussed below. 

Trigger particles are defined as all charged particles originating from the primary vertex, with 
< 2.4 and in a specified p range. The number of trigger particles in the event is denoted 
by Ntrig, which can be more than one per event. Hadron pairs are formed by associating with 
every trigger particle the remaining charged particles with | // 1 < 2.4 and in a specified pj^^°'^ 
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range. The per-trigger-particle associated yield distribution is then defined by: 

' -B(0.0)x4^. (1) 



NtngdAf]dA(p ^ ' ' B{A}],A(p)' 

where At] and are the differences in rj and (p of the pair, respectively. The signal distribution, 
S{Arj, A(p), is the measured per-trigger-particle distribution of same-event pairs, i.e.. 

The mixed-event background distribution, 

B(An,A(p) = -— , (3) 

is constructed by pairing the trigger particles in each event with the associated particles from 10 
different random events, excluding the original event. The symbol N™''' denotes the number of 
pairs taken from the mixed event. The background distribution is used to account for random 
combinatorial background and pair-acceptance effects. The normalization factor B(0, 0) is the 
value of B ( A?/, A(p) at Af/ = and A(p = (with a bin width of 0.3 in A?/ and ti/ 16 in Ai^), repre- 
senting the mixed-event associated yield for both particles of the pair going in approximately 
the same direction, thus having full pair acceptance. Therefore, the ratio B{0,0) / B{Af], A(p) 
is the pair-acceptance correction factor used to derive the corrected per-trigger-particle associ- 
ated yield distribution. Equation ^ is calculated in 2 cm wide bins of the vertex position (Zvtx) 
along the beam direction and averaged over the range |zvtx| < 15 cm. To maximize the statis- 
tical precision, the absolute values of Arj and Acp are used to fill one quadrant of the {Ari, Acp) 
histograms, with the other three quadrants filled (only for illustration purposes) by reflection. 
Therefore, the resulting distributions are symmetric about {Arj,A(p) = (0,0) by construction. 

Each reconstructed track is weighted by the inverse of the efficiency factor, etrkl'^/Pi)/ as a 
function of the track's pseudorapidity and transverse momentum. The efficiency weighting 
factor accounts for the detector acceptance A{r], pj), the reconstruction efficiency E{t], pj), and 
the fraction of misidentified tracks, F{t], pj), 

AE 

^trkiV'Pr) = F- (4) 



Studies with simulated Monte Carlo (MC) events show that the combined geometrical accep- 
tance and reconstruction efficiency for the primary-track reconstruction reaches about 60% for 
the 0-5% most central PbPb collisions at pj > 2 GeV/ c over the full CMS tracker acceptance 
(I //I < 2.4) and 65% for \f]\ < 1.0. The fraction of misidentified tracks is about 1-2% for 
\r]\ < 1.0, but increases to 10% at |?/| ~ 2.4. The weighting changes the overall scale but 
not the shape of the associated yield distribution, which depends on the ratio of the signal to 
background distributions. 

A closure test of the track-weighting procedure is performed on HYDJET Ii33j (version 1.6) 
MC events. The efficiency-weighted associated yield distribution from reconstructed tracks is 
found to agree with the generator-level correlation function to within 3.3%. In addition, sys- 
tematic checks of the tracking efficiency, in which simulated MC tracks are embedded into 
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3 Results 



data events, give results consistent with pure HYDJET simulations to within 5.0%. The track- 
ing efficiency also depends on the vertex z position of the event. However, this dependence 
is negligible in this analysis, and its effects are taken into account in the systematic uncer- 
tainty by comparing the efficiency-corrected correlation functions for two different Zvtx ranges, 
|Zvtx I < 15 cm and |zvtx | < 5 cm, which are found to differ by less than 2.2%. Additional uncer- 
tainties due to track quality cuts are examined by loosening or tightening the track selections 
described previously, and the final results are found to be insensitive to the selections to within 
2.0%. An independent analysis, using a somewhat different but well-established methodol- 
ogy | |34ll35| in constructing the mixed-event background is performed as a cross-check, where 
10 trigger particles from different events are selected first and combined to form a single event, 
and then correlated with particles from another event. It yields results within 2.9-3.6% of the 
default values (3.6% for |A?7| < 1 and 2.9% for 2 < |A?^| < 4), with a slight dependence on Af/ 
and A^. The other four sources of systematic uncertainty are largely independent of ^rj and 

^<p. 

Table [l] summarizes the different systematic sources, whose corresponding uncertainties are 
added in quadrature becoming the quoted systematic uncertainties of the per-trigger-particle 
associated yield. 



Table 1: Summary of systematic uncertainties. 



Source 


Systematic uncertainty of 
the per-trigger-particle associated yield (%) 


Tracking weighting closure test 


3.3 


Tracking efficiency 


5.0 


Vertex dependence 


2.2 


Track selection dependence 


2.0 


Construction of the mixed-event background 


2.9-3.6 


Total 


7.3-7.6 



3 Results 

The measured per-trigger-particle associated yield distribution of charged hadrons as a func- 
tion of I A?/ 1 and \^(p\ in the 0-5% most central PbPb collisions at = 2.76 TeV is shown in 

Fig. [l^ for trigger particles with 4 < p^'^ < 6GeV/c and associated particles with 2 < p|**°'^ < 

4 GeV/ c. To understand the effects of the hot, dense medium produced in the collisions, this 
distribution can be compared to that from a pythiaS MC simulation Il36l (version 8.135) of pp 
collisions at a/s = 2.76 TeV, shown in Fig. This transverse momentum range, one of several 
studied later in this paper, is chosen for this figure since it illustrates the differences between 
correlations from PbPb data and PYTHIA8 pp MC events most clearly. The main features of the 
simulated pp MC distribution are a narrow jet-fragmentation peak at (Af/, A^) (0,0) and 
a back-to-back jet structure at |A^| = n, but extended in A//. In the 0-5% most central PbPb 
collisions, particle correlations are significantly modified, as shown in Fig. [l^. The away-side 
pairs (A^ n) exhibit a correlation with similar amplitude compared to PYTHIAS, although 
the structure in PbPb data is much broader in both Ai/) and Af/ so that it appears almost flat. 
On the near side (A(^ « 0), besides the common jet-like particle production in both pp and 
PbPb at (A)/, A^) ^ (0, 0) due to jet fragmentation, a clear and significant ridge-like structure 
is observed in PbPb at A^ ?a 0, which extends all the way to the limit of the measurement of 
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In relativistic heavy-ion collisions, the anisotropic hydrodynamic expansion of the produced 
medium is one possible source of long-range azimuthal correlations, driven by the event-by- 
event initial anisotropy of the collision zone [7, 37J. For non-central collisions, these correlations 
are dominated by the second-order Fourier component of the |A(^| distribution, usually called 
elliptic flow or V2. Measurements of dihadron correlations at RHIC have frequently attempted 
to subtract or factorize the elliptic flow contribution based on direct V2 measurements, in order 
to reveal other features of particle correlations that may provide insight into the interactions 
between the jets and the medium. However, recent theoretical developments indicate that 
the interplay between initial-state fluctuations and the subsequent hydrodynamic expansion 
gives rise to additional Fourier components in the azimuthal particle correlations |25l l26l l38l - 
1^5. These components need to be treated on equal footing with the elliptic flow component. 
In particular for the 0-5% most central PbPb collisions, the elliptic flow contribution to the 
azimuthal correlations is not expected to be dominant I43|. Therefore, the original unsubtracted 
correlation functions are presented in this paper, containing the full information necessary for 
the comparison with theoretical calculations. 




Figure 1: Two-dimensional (2-D) per- trigger-particle associated yield of charged hadrons as a 
function of |A?/| and |A^| for 4 < p^"^ < 6GeV/c and 2 < pf^^" < 4GeV/c from (a) 0-5% most 
central PbPb collisions at ./s^ = 2.76 TeV, and (b) PYTHIA8 pp MC simulation at = 2.76 TeV. 



3.1 Associated Yield Distributions versus 

To quantitatively examine the features of short-range and long-range azimuthal correlations, 
one dimensional (1-D) correlation functions are calculated by averaging the 2-D distribu- 
tions over a limited region in Arj from A^^min to Arjmax- 

-dArj. (5) 



Ntrig dA(p A//max " Aj/min J Aij^in Mrig dAt]dA(p 

The results of extracting the 1-D A^ correlations for the 0-5% most central PbPb collisions are 
shown in Figs. |2] and |3] The associated yield per trigger particle in the range of 2 < pj^^°^ < 

4GeV/c is extracted for five different pj'^ intervals (2-4, 4-6, 6-8, 8-10, and 10-12 GeV/c) and 
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3 Results 



two ranges in Arj. Figure |2] gives the short-range pseudorapidity result, i.e., averaged over the 
region \ At]\ < 1. Figure [sjshows the same comparison for long-range correlations, i.e., averaged 
over the region 2 < |Aj/| < 4. A comparison to PYTHIA8 pp MC events at = 2.76 TeV is 
also shown, with a constant added to match the PbPb results at = 1 in order to facilitate 
the comparison. In this projection, only the range < Ai^ < tt is shown, as the A^ correlation 
function is symmetric around A(p = by construction. 



< |Ati| < 1 



2 < p^'""" < 4 GeV/c CMS J L dt = 3.1 nb"' 



CM 

"D 



^ 6.5 



2<p/<4GeV/c 

PbPb^s^ = 2.76 TeV, 0-5% 
PYTHIAS pp>is = 2.76T6V 
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Figure 2: Short-range (|A//| < 1) per-trigger-particle associated yields of charged hadrons as 
a function of |A^| from the 0-5% most central PbPb collisions at = 2.76 TeV, requiring 

2 < pf^°'^ < 4 GeV/c, for five different intervals of pj'^. The pythiaS pp MC results (solid 
histograms) are also shown, shifted up by a constant value to match the PbPb data at A^ = 1 
for ease of comparison. The error bars are statistical only and are too small to be visible in most 
of the panels. The systematic uncertainty of 7.6% for all data points is not shown in the plots. 



2 < |Ati| < 4 



2 < p^''"" < 4 GeV/c CMS I L dt = 3.1 nb-' 




m M m M m 

Figure 3: Long-range (2 < \Arj\ < 4) per-trigger-particle associated yields of charged hadrons 
as a function of |A^| under the same conditions as in Fig. [2] The systematic uncertainty of 7.3% 
for all data points is not shown in the plots. Dashed lines show the fits by the first five terms in 



the Fourier series, as discussed in Section 3.3 



The panels of Figs. [2] and [s] show a modification of the away-side associated yield for \A(p\ > 
n/ 2 in central PbPb collisions that is not present in the PYTHIA8 pp MC simulation. This modi- 
fication is most pronounced for lower p values and is characterized by a similarly broad dis- 
tribution as seen in lower-energy RHIC measurements, where it was attributed to jet quenching 
and related phenomena L3}i5J. The near-side associated yield in PbPb data includes the con- 
tributions from both the jet-like peak and the ridge structure seen in Fig. [l^. Thus, it is not 
directly comparable to pythiaS pp MC events, where the ridge component is absent. A more 
quantitative comparison of the jet-like component between PbPb data and pythiaS pp MC 
events will be discussed below. 

For azimuthal correlations at large values of Arj (Fig. |3|, a clear maximum at A^ is ob- 
served, which corresponds to the ridge structure seen in the 2-D distribution of Fig. [T^. This 
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feature of the long-range azimuthal correlation function is not present in the PYTHIAS pp MC 
simulation for any p bin. In PbPb, the height of the ridge structure decreases as pj^^ increases 
(Fig. ij and has largely vanished for pj^^ 10-12 GeV/c. The diminishing height of the ridge 

with increasing pj^^ was not evident from previous measurements at RHIC in AuAu collisions 
presumably because of lack of events at high pj. 

3.2 Integrated Associated Yield 

The strengths of the jet peak and ridge on the near side, as well as their dependences on Arj and 
Pj^^, can be quantified by the integrated associated yields. In the presence of multiple sources 
of correlations, the correlation of interest is commonly estimated using an implementation of 
the zero-yield-at-minimum (ZYAM) method [17J. However, as mentioned previously, the pos- 
sible contribution of elliptic flow is not taken into account as it is not the dominant effect for the 
most central PbPb collisions considered here. The ZYAM method is implemented as follows. 
A second-order polynomial is fitted to the |A(^| distributions in the region 0.5 < |A^| < 1.5. 
The location of the minimum of the polynomial in this region is denoted as A^zyam- Using 
the position of the minimum, the associated yield is then calculated as the integral of the \A(p\ 
distribution minus its value at A(/)zyam between | A^| = and A^zyam- The uncertainty on the 
minimum level obtained by the ZYAM procedure, combined with the uncertainty arising from 
the choice of fit range in |A^|, results in an uncertainty on the absolute associated yield that is 
constant with a value of 0.012 over all Arj and pj^^ bins. 
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Figure 4: Integrated near-side (|A^| < A^zyam) associated yield for 4 < pj^ < 6GeV/c and 
2 < p|^*°'^ < 4 GeV/ c, above the minimum level found by the ZYAM procedure, as a function of 
\Ar] I for the 0-5% most central PbPb collisions at ^/s^ = 2.76 TeV. The error bars correspond to 
statistical uncertainties, while the brackets denote the systematic uncertainties. The solid line 
shows the prediction from the PYTHIAS simulation of pp collisions at i/s = 2.76 TeV. 



The ZYAM procedure enables the direct extraction of integrated yields of the A</)-projected 
distributions in a well-defined manner. Figure |4] shows the resulting near-side associated yield 
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3 Results 
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Figure 5: Integrated near-side (|A^| < AcpzYAu) associated yield above the minimum level 
found by the ZYAM procedure for (a) the short-range jet region (|A?/| < 1) and (b) the long- 
range ridge region (2 < \Arj \ < A), as a function of Pj'^ in the 0-5% most central PbPb collisions 
at ^/s^ = 2.76 TeV with 2 < pf^°'^ < 4 GeV/c. The error bars correspond to statistical uncertain- 
ties, while the brackets around the data points denote the systematic uncertainties. The solid 
line shows the prediction from the pythiaS simulation of pp collisions at = 2.76 TeV. 



as a function of | A// 1 (in slices of 0.6 units) in both the central PbPb collisions and PYTHIAS pp 
MC. The near-side associated yield for pythiaS pp shows a strong peak at \At]\ = 0, which 
corresponds to the expected correlations within jets. This near-side peak diminishes rapidly 
with increasing |A//|. The PbPb data also exhibit a jet-like correlation peak in the yield for 
small \Ai]\, but in contrast, the PbPb data clearly show that the ridge extends to the highest 
\Arj \ values measured. 

Figure |5] presents the integrated associated yield for the jet region (|A?/| < 1) and the ridge 
region (2 < |Aj/| < 4) with 2 < pf^°^ < 4GeV/c, as a function of p^'^ in the 0-5% most 
central PbPb collisions. The ridge-region yield is defined as the integral of the near side in 
long-range A(p azimuthal correlation functions (Fig. |3|, while the jet-region yield is determined 
by the difference between the short- and long-range near-side integral, as the ridge is found 
to be approximately constant in Arj (Fig. |4j|. While the jet-region yield shows an increase with 
Pj'^ due to the increasing jet transverse energy, the ridge-region yield is most prominent for 

2 < pj'^ < 6GeV/c and tends to drop to almost zero when pj^^ reaches 10-12 GeV/c. The 
bars in Fig. |5] correspond to the statistical uncertainties, while the brackets around the data 
points denote the systematic uncertainties, which are dominated by the tracking performance, 
as discussed earlier. Results from the pythiaS pp MC simulation, displayed as solid lines in 
Fig. 5 are consistent with zero for all Pj'^ bins in the ridge region and have qualitatively the 

same trend with pj^^ in the jet region as the data. It has been seen from previous studies that 
the ridge is absent in both minimum bias pp data and pythiaS MC simulations. However, the 
PYTHIAS program does not fully describe the jet-like correlations observed in pp collisions at 
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= 7TeV pTi . A quantitative comparison of the correlations in PbPb and pp collisions will 
therefore be deferred until pp results at -^s = 2.76 TeV become available. 
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Figure 6: Fourier coefficients, Vf, V3, y|, and extracted as functions of p^'^ for 2 < 
passoc ^ 4GeV/c for the 0-5% most central PbPb collisions at ^Jj:;^ = 2.76 TeV. The error bars 
represent statistical uncertainties only. The solid lines show the predictions from the PYTHIAS 
simulation of pp collisions at ^/s = 2.76 TeV. 



3.3 Fourier Decomposition of Distributions 

Motivated by recent theoretical developments in understanding the long-range ridge effect 
in the context of higher-order hydrodynamic flow induced by initial geometric fluctuations, 
such as the "triangular flow" effect Il26l l38-42l, an alternative way of quantifying the observed 
long-range correlations is investigated in this paper. The 1-D A^-projected distribution for 
2 < \Ar] I < 4 is decomposed in a Fourier series as 



1 dNP^i'' N, 



Ntrig dA(p 27r 1^ — ^ 

where Nassoc represents the total number of hadron pairs per trigger particle for a given | Arj \ 
range and (pj'^, bin. The 1-D projections displayed in Fig. [i are fitted by the first 

five terms in the Fourier series, the resulting fits being shown as the dashed lines in Fig. |3] 
The data are well described by the fits. Figure |6] presents the first five Fourier coefficients 
from the fit as functions of p^'^ for 2 < pf^°'^ < 4GeV/c for the 0-5% most central PbPb 
collisions and for the PYTHIA8 pp MC simulation. The PYTHIA8 results are scaled by the ratio 
of Nlg™^^/N^^J^ in order to remove the trivial effect of the multiplicity dependence. The error 
bars are statistical only, while the systematic uncertainties are found to be negligible because 
the Fourier coefficients characterize the overall shape of the correlation functions, and thus are 
not sensitive to the absolute scale. The coefficients from the fit are also found to be largely 
independent of each other (correlation coefficients typical below 5%). 

All the Fourier coefficients found from fitting the PbPb data show a similar dependence on 
Pj'^ except for the Vf term, which contains an additional negative contribution that grows 

toward higher Pj'^. This negative component is consistent with a contribution from mo- 
mentum conservation or back-to-back dijets [44]. If the observed correlation is purely driven 
by the single-particle azimuthal anisotropy arising from the hydrodynamic expansion of the 
medium |45|, the extracted V,{ components would be related to the flow coefficients v„ (i.e., 
V2 for anisotropic elliptic flow) via ~ vf,^^ x f where vj^^ and v^^'^ are the flow co- 
efficients for the trigger and associated particles p6|. The flow coefficients, and especially 



l+X]2yfcos(nA</))l, (6) 
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the higher-order terms, are sensitive to the initial conditions and viscosity of the hot, dense 
medium l|38l l46ll . This Fourier analysis serves as an alternative way of quantifying the az- 
imuthal correlation functions in heavy-ion collisions and potentially provides new insights in 
understanding the origin of the ridge effect in these collisions. Analysis of larger data samples 
in terms of p^"^ and p|®®°'^, as well as centrality, would allow a more detailed comparison to 
theoretical calculations of hydrod}mamics, for which the Fourier components provide a con- 
cise description of the data. 

4 Summary 

The CMS detector at the LHC has been used to measure angular correlations between charged 
particles in At] and Acp up to | At] \ 4 and over the full range of in the 0-5% most central 
PbPb collisions at a/s^ = 2.76 TeV. This is the first study of long-range azimuthal correlations 
over a large difference in pseudorapidity in PbPb interactions at the LHC energy. The extracted 
2-D associated yield distributions show a variety of characteristic features in heavy-ion colli- 
sions that are not present in minimum bias pp interactions. Short- and long-range azimuthal 
correlations have been studied as a function of the transverse momentum of the trigger parti- 
cles. The observed long-range ridge-like structure for approximately equal azimuthal angles 
(A^ 0) is most evident in the intermediate transverse momentum range, 2 < pj^^ < 6 GeV/ c, 

and decreases to almost zero for pj'^ above 10-12 GeV/c. A qualitatively similar dependence 
of the ridge on transverse momentum has also been observed in high-multiplicity pp events, 
indicating a potentially similar physical origin of the effect. A Fourier decomposition of the 
1-D A(^-projected correlation functions in the ridge region (2 < \Af]\ < 4) has been presented. 
This alternative way of quantifying the correlation data provides valuable information to test a 
wide range of theoretical models, including recent hydrodynamic calculations of higher-order 
Fourier components. The very broad solid-angle coverage of the CMS detector and the statisti- 
cal accuracy of the sample analyzed in this paper provide significantly improved observations 
of short- and long-range particle correlations over previously available measurements. 
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